Ground and excited states of spinor Fermi gases in tight waveguides and the 
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The ground and excited states of a one-dimensional (ID) spin-| Fermi gas (SFG) with both attrac- 
tive zero-range odd-wave interactions and repulsive zero-range even-wave interactions are mapped 
exactly to a ID Lieb-Liniger-Heisenberg (LLH) model with delta-function repulsions depending on 
isotropic Heisenberg spin-spin interactions, such that the complete SFG and LLH energy spectra 
are identical. The ground state in the ferromagnetic phase is given exactly by the Lieb-Liniger (LL) 
Bethe ansatz, and that in the antiferromagnetic phase by a variational method combining Bethe 
ansatz solutions of the LL and ID Heisenberg models. There are excitation branches corresponding 
to LL type I and II phonons and spin waves, the latter behaving quadratically for small wave number 
in the ferromagnetic phase and linearly in the antiferromagnetic phase. 
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If an ultracold atomic vapor is confined in a de Broglie 
wave guide with transverse trapping so tight and tem- 
perature so low that the transverse vibrational excita- 
tion quantum Hu> is larger than available longitudinal 
zero point and thermal energies, the effective dynam- 
ics becomes one-dimensional (ID) [1, 2], a regime cur- 
rently under intense experimental investigation [3, 4]. 
Confinement-induced ID Feshbach resonances (CIRs) 
reachable by tuning the ID coupling constant via 3D 
Feshbach scattering resonances [5] occur for both Bose 
[1, 6] and spin- aligned Fermi gases [7], allowing tuning of 
the effective ID coupling constants to any desired values. 
Near the CIRs the ground states have strong short-range 
correlations not representable by effective field theories, 
and such systems have recently become the subject of 
extensive theoretical and experimental investigations be- 
cause of their novelty and the recent experimental realiza- 
tion [4] of the ID gas of impenetrable point bosons, now 
known as the Tonks-Girardeau (TG) gas. The Fermi- 
Bose (FB) mapping method, introduced in 1960 [8, 9] 
to solve the TG problem exactly, and recently [7, 10-13] 
extended to a mapping between a ID Bose gas with delta- 
function interactions of any strength [Lieb-Liniger (LL) 
gas [14]] and a spin-aligned Fermi gas, will be generalized 
here to the Fermi gas with unconstrained spins. 

The exact ground state of the "fermionic TG" (FTG) 
gas, a spin- aligned ID Fermi gas with infinitely strongly 
attractive zero-range odd-wave interactions, has been de- 
termined recently [7, 11-13] by FB mapping to the ideal 
Bose gas, and its highly-correlated ground state was 
shown to possess superconductive off-diagonal long-range 
order (ODLRO) [15, 16]; this system models a magnet- 
ically trapped ultracold gas of spin-i fermionic atoms 
with very strong odd-wave interactions induced by a p- 
wave Feshbach resonance. If optically trapped instead, 
such a spinor Fermi gas (SFG) has richer properties, 
since its spins are unconstrained. In such a case, with 
both spatial even-wave interactions associated with spin 



singlet scattering and odd-wave interactions associated 
with spin triplet scattering, one can vary the ratio of 
the two coupling constants by Feshbach resonance tun- 
ing of the odd-wave one, leading to a rich phase diagram 
of ground state spin, which we have determined recently 
[17]. However, exact ground and excited states in the 
antiferromagnetic region of the phase diagram are un- 
known except for limiting cases where one of these in- 
teractions is negligible. The purpose of this paper is to 
study the more general problem by a generalized space- 
spin FB mapping to a ID Lieb-Liniger-Heisenberg (LLH) 
model of particles with delta-function repulsions depend- 
ing on Heisenberg spin-spin interactions, followed by a 
variational treatment combining Bethe ansatz solutions 
of the LL and ID Heisenberg models. 

Even and odd-wave SFG interactions: The SFG is an 
ultracold vapor of fermionic atoms in ID in two different 
hyperfine states, representable as spin-i fermions. One 
can represent the TV-atom states by space-spin wave func- 
tions iPf{xi, ft; • • • ',%n, &n) where the ID spatial 
coordinates and the spin z-component variables take on 
values ±| labelled by f and J,. This system is assumed to 
be optically trapped on a ring of circumference L, so the 
Xj are measured circumferentially and periodic boundary 
conditions with periodicity length L are imposed. The 
Hamiltonian with LL units h = 2m = 1 in a limit of 
zero-range interactions is [11-13, 17] 
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Here Xjt — Xj — xg, d is the Dirac delta function [14], 

( \t> * s P r °j ector onto tnc subspaces of singlet 
(triplet) functions of the spin arguments (tjj, <jg) for fixed 
values of all other arguments, and v° D is a strong, at- 
tractive, odd-wave ID interaction of range xq — > 0+ 
(1D analog of 3D p-wave interactions) occuring near the 
CIR [7, 11—13]. It is the zero-range limit of a deep 
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and narrow square well of depth Vq and width 2xq with 
xq — ► 0+ and Vo — > oo such that Vox^ approaches a finite 
limit, implying an attractive interaction stronger than a 
delta function, which induces a wave-function disconti- 
nuity in the zero-range limit [10]. Taking the limit such 
that = (tt/2x )[1 + (2 / V) 2 [ (x / 'a° 1D )] generates a 

two-body wave function satisfying the contact condition 
ipF(xje = 0+) = -ip F (xje = 0-) = -a° D ip F (xje = 0±) 
where a° D < is the ID odd- wave scattering length and 
the prime denotes differentiation [7, 11-13]. The even- 
wave ID coupling constant gf D in (1) is related to the 
even-wave scattering length a\ D < derived [1] from 
3D s-wave scattering by g\ D = —A/a\ D and the even- 
wave contact condition is the usual LL one [14] ip' F (xjt = 
0+) = -ip' F (xje = 0-) = -{aloT^Fixjt = 0±). 

The (j, €) singlet subspace is spanned by the bispinor 
wave function X s { <J j, = ^(hli - Ijh) and the 
triplet subspace by x\ =TjT<, Xo = 7f ( L?-U + Ijh), and 
xL\ =ljle- Since tpp is antisymmetric under combined 
space-spin exchanges (xj,<jj) <-» (a^,0>), the singlet is 
the part of ipF which is even in ir^ as — > and the 
triplet is the part which is odd. In [11-13, 17] we ac- 
counted for this entirely by definition of the even and 
odd-wave pseudopotentials derived from 3D s-wave and 
p-wave scattering, and omitted the singlet and triplet 
projectors. However, here we study the SFG by a map- 
ping which converts the odd waves into even ones while 
leaving the even ones unchanged, so it is essential to ad- 
join these projectors so as not to lose track of this as- 
sociation. They are related to the isotropic Hciscnberg 
spin-spin interaction Sj ■ Si by P- f = \ — Sj ■ Sg and 

jt = I J ' e wnere are the usual spin vector 
operators. 

SFG-LLH mapping: Assume that the system is con- 
tained in a ID ring trap of circumference L, with peri- 
odic boundary conditions of periodicity length L. The 
spinor Fermi gas wave functions and Hamiltonian can 
be mapped to those of a Lieb-Liniger-Heisenberg model, 
with wave functions iPllh{xi, <7i; • • • ; xn, <tn) invariant 
under cyclic permutations [x\, o\\ X2, <Ji\ • • ■ ; xn 1 o n) — > 
(x2,cr 2 ; - ■ ■ ;xN,cr N ]Xi,ai) around the ring, by multiply- 
ing the part of the SFG wave function ipp which is odd 
in Xji as Xjg — > by e(xjg) where e(x) = +1 (—1) for 
x > (x < 0), and e(0) = 0, while leaving the even 
part unchanged; this generalizes the original FB map- 
ping of [8]. The iV-particle configuration space decom- 
poses into TV! permutation sectors, each defined by one 
particular ordering of numerical values of x\, ■ ■ ■ , xn- In 
analogy with the LL model [14] it is sufficient to define 
our spinor fermion wave functions tpp only in the funda- 
mental permutation sector 9ti wherein < x\ < x 2 < 
■ ■ ■ < xn < N and to specify boundary conditions on the 
boundary of 9ti, which is an (N — l)-dimensional hyper- 
surface Si. This boundary can be crossed in each of N 
orthogonal directions, each defined by a particular coor- 



dinate Xj crossing its right-hand neighbor Xj + \ so that 
one enters a neighboring permutation sector wherein 
xi < x 2 < ■ ■ ■ < Xj-i < x j+ i < x,j < Xj +2 < ■ ■ ■ < xn- 
One obtains a well-posed problem for determination of 
the energy eigenfunctions ipFa and eigenvalues E a by 
requiring that tpFa satisfy the free-particle Schrodingcr 
equation in the interior of 9li plus boundary conditions 
relating its value and normal gradient on Si; ipFa can 
be uniquely extended to all N\ permutation sectors by 
space-spin antisymmetry. 



Since the spatially odd (even) part of ipFa is associ- 
ated with the spin triplet (singlet) subspace, the desired 
mapping tpp iPllh = UipF can be implemented by 
a unitary mapping operator U = UnUn-i • • • U± where 
each factor for 1 < j < N — 1 is Uj = Pjj+i + e ( x j+i — 
Xj)Pj j +1 . Uj is unitary and self- inverse and changes the 
component of i\)p odd in Xjj + \ around Xjj + \ = into an 
even one while leaving the even component unchanged. 
Un is a special case required to close the ring periodi- 
cally: Un = Pjtf i + e(xi + L — xn)P n i- Since Xj < Xj + \ 
in *Hi, the e(xj+i — Xj) factors are all unity there, and 
in £fti the mapping reduces to the identity. The term 
e(xj + i — xj) in Uj is to be interpreted as a prescription 
for extension to the adjacent permutation sector JHi :J . In 
91 1 j all factors Ug with £ ^ j are unity, verifying that U 
maps i/if to a LLH wave function ipLLH differing from 
ipF only by extension from EHi into adjacent sectors £Hi,j 
by multiplication of the component of ipF odd in Xjj+i 
around Xjj+i = by e(xj + \ — Xj). 



When extended to all permutation sectors ipLLH is nei- 
ther bosonic nor fermionic under combined space-spin ex- 
changes (xj, Uj) <-> (xi, ag); it has both Bose-like compo- 
nents which are space-even and spin-even, and Fermi-like 
components which are space-even and spin-odd, but no 
space-odd components. The mapping converts the odd- 
spatial-wave, spin triplet contact discontinuities of ipF 
into LL contact cusps where ipp is continuous and its 
derivative changes sign, while retaining the original LL 
even-spatial- wave, spin singlet cusps. As in [11-13] this 
changes the odd spatial wave, spin triplet contact condi- 
tion tpF{xje = 0+) = -ip F (xj£ = 0-) = -a° lD il) F (xjg = 
0±) of the SFG into an even spatial wave, spin triplet 
contact condition tp' LLH (xj£ = 0+) = —^P'llh{ x 3^ = 
0— ) = —{a c l D )~ 1 il)[ J [ J H{xji = 0±) for the LLH gas with 
the same scattering length a° D . Then both singlet and 
triplet contact conditions are of the usual form generated 
by LL delta functions [14] and it follows that ipLLHa is 
an eigenstate of Hllh with energy E a if ipp a is an eigen- 
state of Hp with energy E a and conversely, where after 
substituting P? e = \ — Sj ■ Sg and P* e — § + Sj • one 
finds 
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V 4- d 

H LLH = -J2 d ^+ E [-^-+2{c-(S)S j * e ]5(x je )] 

3=1 l<j<£<N 

(2) 

with coupling constants given in terms of the odd and 
even-wave SFG scattering lengths a° D and a\ D by c = 
2/K D | and c' = = 2/\a\ D \. 

Ground state phase diagram: Assume that N is even. 
We previously determined [17] the phase diagram of total 
ground-state spin 5 in the (c, d) plane, finding (a) for 
c < d the ground state is ferromagnetic with S = N/2; 
(b) for c > d it is antifcrromagnetic with 5 = 0; (c) 
for c = d it is highly degenerate, all (integer) values of 
5 from to | and 5 Z from -| to | having the same 
energy, (c) now follows trivially from Eq. (2); if c = d the 
Sj • terms in (2) vanish, so the entire energy spectrum 
is spin- independent. 

Ferromagnetic phase: In this phase, which occurs for 
c < d, the ground state is totally spin-aligned (5 = N/2), 
there are no quantum fluctuations of Sj ■ Sg which may 
therefore be replaced by its expectation value j, and 
the square bracket in (2) reduces simply to 2c5(xji), the 
usual LL interaction, independent of d . Then the LL 
Bcthe ansatz [14] yields the exact ground state for all c, 
as in [7, 11-13]. As c increases from zero for fixed d, the 
system remains in the ferromagnetic phase until c reaches 
d, where it undergoes a quantum phase transition to the 
5 = antiferromagnetic phase. This is illustrated in 
Fig. 1, where the solid line shows the ground state en- 
ergy Eq(c) in the ferromagnetic phase [14, 18] and the 
dashed lines for fixed values of d show the energy in the 
antiferromagnetic phase. 

Antiferromagnetic phase: Now assume c > d , in which 
case the ground state is highly nontrivial, has 5 = 0, 
and depends on both c and d . As before, it is sufficient 
to restrict ourselves to the sector 9ti: < x\ < X2 < 
■ ■ ■ < xn < L. One can obtain an upper bound on the 
ground state energy by taking a variational state which 
is a product of a purely spatial state V'space^i, ■ ■ ■ ,xn) 
by a purely spin state V'spin(c r i, • • • , ctjv). In spite of 
its product structure, this trial state incorporates some 
space-spin correlation implicitly via the projectors P? e = 

j—Sj- S e and P? e = | + Sj ■ Sg. Taking Vspin to be the 
ground state of the isotropic, antiferromagnetic Heisen- 
berg Hamiltonian Hueis — ^2^=1 Sj ' Sj+i + Sjv ■ Si, for 
which in the thermodynamic limit -Eo.Heis/A" = j — In 2 
[21], reduces the Hamiltonian to an LL model [14, 22], 
yielding an upper bound Eq(c,c') < -Bo.ll(cll) with 
cll = c(l — In 2) + c'ln2. The product state Vs paceV'spin 
is exact in the ferromagnetic phase (where Sj -Sj+i = \), 
on the phase boundary c = d, and when c = [23, 24], 
so it may be a good approximation throughout the (c, d) 
plane. This variational ground state energy, obtained 
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FIG. 1: Ground state energy of the ID spinor Fermi gas as a 
function of c = 2/\a° D \ for fixed values of c' = 2/|af D |, where 
a° D and a\ D are odd and even-wave ID scattering lengths. 
Solid line: Exact ground state energy in the ferromagnetic 
phase, stable only for c < c'. Dashed lines: Variational 
ground state energy in the antiferromagnetic phase, stable 
only for c > c', for c'/n = 0, 2, 5, and 10 (bottom to top). 

from the LL ground-state energy [18], is shown by the 
dashed lines in Fig. 1. 

Excitations: Both phonons (number density waves) 
and spin waves (spin density waves) occur in this system. 
Since a space-spin product variational approximation was 
used for the ground state, the energies of these excitations 
will be evaluated using variational states V'spaceV'spin, but 
now with either V'spacc or i/'spin excited. For one-phonon 
excitations i/'spacc is an exact one-phonon LL excited 
state [25] and V'spin is the exact ground state of the ID 
Heisenberg model. Consider first such excitations from 
the ferromagnetic phase (c < d). There are two such 
modes: Hole excitations (Lieb's type II excitations [25]) 
correspond to promotion of a particle from quasimomcn- 
tum K — k inside the LL-Fermi sea (Fig. 2 of LL [14]) 
to just above the top (quasimomcntum K) of the sea, 
whereas umklapp excitations (Lieb's type I excitations) 
correspond to promotion of a particle from —K to K + k. 
Both have the same initial slope corresponding to a sound 
speed given, assuming LL units h = 2m = 1, by Lieb's 

Eqs. (1.4) and (1.5), v s = 2^J /j, — Here fi is the 

LL ground state chemical potential /i = n 2 (3e — ■ye), the 
dot denotes differentiation with respect to 7, n is the ID 
atom number density, ef/y) = n~ 2 (E /N), E is the LL 
ground state energy, and 7 = c/n where c is the coupling 
constant in the LL interaction 2c5(xje) [25]. The type I 
and II curves as a function of momentum are shown in 
Lieb's Figs. 1-3. The corresponding phonon excitations 
from the antiferromagnetic phase (c > d) differ only by 
replacement of c by cll = c(l — In 2) + d In 2. 

Determination of spin wave energies is more delicate. 
The total energy is the expectation value of (2) in a 
space-spin product state V'spaceV'spin in which V'space is the 
ground state and V'spin is a one spin wave excited state. 
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FIG. 2: Solid line: ef orro (g)/(c' — c)n A e'(c/n). Dashed line: 

eantiferro(g)/ (c - d )n A d (c LL / n) . 

The spin wave energy enters only indirectly via the ex- 
pectation values (Sj ■ Sj + i) which affect the LL coupling 
constant cll — Zc \ c + (c — d){Sj ■ Sj+i). Consider first 
spin- wave excitations in the ferromagnetic phase (c < d), 
for which {Sj ■ 5j+i)o = \ for all j, and Cll — c. In the 
excited state it is still independent of j, but now {Sj • 
Sj+i) = {Sj ■ Sj+i)o + A(q) where A(q) is the change due 
to the presence of a spin wave of wave number q, given 
by A(q) = -en cis (q)/N [26] with e He is(?) the Heisenberg 
spin wave energy. For Hn cls = — . Sj ■ Sj+i one has 
[27] e H cisO?) = 2(1- cos q) and A{q) = -2A^- 1 (1 -cosq), 
yielding a change —2N~ 1 (c — c')(l — cosg) in the effec- 
tive LL coupling constant from the expectation value 
of Eq. (2). Then Taylor expansion about the ground 
state energy of Hllh yields the SFG spin wave energy 
£fcrro(<z) = 2(c' — c)(l — cosq)n 3 e'(j) where e'(-f) is the 
derivative of e(-f) = E /n 2 N, 7 = c/n, and n is the ID 
number density [14, 18]. ef e rro(?) is quadratic for q, as 
is also the case for isospin waves in a ID two-component 
Bose gas [19, 20]. The antiferromagnetic case is almost 
the same, except that the change in {Sj ■ Sj+i) from 
its ground state value ^ — In 2 is now positive and de- 
rived from the antiferromagnetic spin wave energy [28] 
£Hcis(g) = f|sing|, yielding a SFG spin wave energy 
e a ntifcrro(<z) = f (c-c')| sin <7|n 3 e'(7 LL ) with 7 Li = c LL /n 
and c LL = c(l — In 2) + d In 2. This is linear for small q. 
These spin wave energies are shown in Fig. 2; the values 
of c and d enter only in scale factors. 

Conclusions: Ground and low excitation energies of 
the ID spinor Fermi gas with both even and odd-wave 
scattering have been determined by an exact mapping to 
a "Lieb-Liniger-Heisenberg" model which was treated by 
a variational method combining Bethe ansatz solutions of 
the LL and Heisenberg models. The richness of ground 
state phases and excitations suggests that experimental 
study of this system could be quite fruitful. 
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